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Reasons have been found for thinking that the minimum diameter of channels of 
a given length to support superconductivity at room temperature through films 
of oxidised atactic polypropylene (OAPP) is considerably larger than found in 
a model for Bose condensation in an array of nanofilaments [D.M. Eagles, Phil. 
Mag. 85, 1931 (2005)] used previously. This model was introduced to interpret 
experimental results dating from 1988 on OAPP. The channels are thought to 
be of larger diameter than believed before because, for an N-S-N system where 
the superconductor consists of an array of single-walled carbon nanotubes, the 
resistance, for good contacts, is Rq/2N, where N is the number of nanotubes 
and Rq = 12.9 kfi [See e.g. M. Ferrier et al., Solid State Commun. 131, 615 
(2004)]. We assume this would be 2Rq/N for a triplet superconductor with all 
spins in the same direction and no orbital degeneracy, which may be the case for 
nanofilaments in OAPP. Hence one may infer a minimum number of filaments 
for a given resistance. In the present model, the E{K) curve for the bosons is 
taken to be of a Bogoliubov form, but with a less steep initial linear term in the 
dispersion at T c than occurs at low T. This form is different from the simple 
linear plus quadratic dispersion, with a steeper initial slope, used in my 2005 
paper. A combination of theory and experimental data has been used to find 
approximate constraints on parameters appearing in the theory. 
Keywords: one-dimensional sytems - high-temperature superconductivity - oxi- 
dised atactic polypropylene -Bose gas - bipolarons - nanofilament arrays 



1. Introduction 

There were claims in a 47-page deposited paper in Russian in 1988 that super- 
conductivity occurs at room temperature in narrow channels through thin films 
of oxidized atactic polypropylene (OAPP)[l , and claims of this in physics jour- 
nals from 1989 [2] . The first indications of possible superconductivity at room 
temperature came from experiments involving the resistance between pressure 
microcontacts on the top of OAPP films and a conducting substrate. Three 
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types of contact points were found, insulating points, points with medium re- 
sistance, and points with low resistance. (See Fig. 1, which is based on Fig. 
17 of Ref. [I], and has also been published in Ref. 3 ]). The resistances of the 
group of low-resistance points lay in the same range for two widely different 
film thicknesses, making it plausible that the resistance was entirely contact 
resistance, with negligible resistance in the channels between the contacts and 
the substrate at these points. 

More convincing evidence of room-temperature superconductivity came in 
1990 and 1991 from experiments which showed (i) that the conductivity of 
the low-resistance channels was several orders of magnitude larger than that of 
Cu[3J[S], (ii) that the high conductivity was destroyed in pulsed measurements 
by non-thermal means when currents above a critical value were applied[5], 
and (iii) there was a negligible electronic contribution to the thermal conduc- 
tivity, implying a violation of the Wiedemann-Franz law by several orders of 
magnitude[3]. There are some questions in the present author's mind about 
how the results of the four-probe measurements of conductivity reported in Ref. 
[3] came about, but the estimate of the conductivity by indirect methods in Ref. 
[5] seems soundly based. This estimate of conductivity made use of the thermal 
conductivity results of Ref. [BJ , and depended on the lack of decomposition of 
the sample after repeated pulsed measurements, which permitted a deduction of 
conductivity greater than 10 ohm cm , and a transition temperature from 
the highly conducting state under pulsed-current conditions of at least 700 K, 
the decomposition temperature of the polymer. Similar very high superconduct- 
ing T c 's have been inferred by Zhao and coworkers EJ [10] from conductivity 
and other measurements on carbon nanotubes, following a suggestion by Tse- 
bro et al.|llj that there might be superconductivity at room temperature in 
carbon nanotubes. Work claiming possible room-temperature superconductiv- 
ity in highly oriented pyrolytic graphite|12[ 113] may be related. Another claim 
of room-temperature superconductivity in a quasi one-dimensional system has 
been made by Djurek et al.[T2] E3 [IB]. In Djurck's work, the main mate- 
rial involved, Ag 5 P 2 06, contains quasi one-dimensional channels, but in recent 
work [15] [IB] it is found that indications of room-temperature superconductivity 
occur only when pellets of Ag 5 P 2 06 are embedded in an insulating matrix. 

Unusual magnetic properties have also been reported in OAPP and other 
polymers. These are: 1. Metamagnetic transitions at magnetic fields of the 
order of 1000-4000 OerstedpH EH EH [20] in films of both OAPP and poly- 
dimcthylsiloxane (PDMS); 2. Large diamagnetism in some films of OAPPjllEl 
|2"D"] : 3. Occasional spontaneous sudden jumps of films of OAPP towards low-field 
regions in inhomogeneous magnetic fields [18]. The metamagnetism is thought 
to be due either to a transition from an antiferromagnetic to a ferromagnetic 
arrangement [Tl [T71 [2"] or to alignment of ferromagnetic domains of ferromagnetic 
(and superconducting) filaments. We take the second viewpoint in this paper. 
The large diamagnetism and spontaneous jumps to low-field regions in inho- 
mogeneous magnetic fields are thought to be associated with superconducting 
channels which join at their ends to form closed loops [U [2] [TH] [20] [21]. In Ref. 
[21] the large diamagnetism observed in two samples of OAPP [20] has been in- 
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Figure 1: Distribution function for resistance R at room temperature from mi- 
croprobes to a metallic substrate: (a) For a sample of oxidised atactic polypropy- 
lene (OAPP) of thickness 50 /im, (b) for OAPP with thickness 0.3 fim, and (c) 
for the metallic substrate without any polymer film. [Reproduced from Fig. 1 
of Rcf.|3 1, which was obtained approximately from Fig. 17 of Ref.pQ] 
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terpreted in detail as due to induced currents in closed loops, using a model with 
loops of two different radii as an approximation to a continuous distribution. 

Almost no attempts to repeat and extend the results mentioned above have 
been reported in the literature. At one point, some members of the Ioffe Institute 
group who reported results supporting superconductivity at room temperature 
in 1990, decided that the effects observed were due to formation of a conduct- 
ing filament of the electrode material in high electric fields [22, 23;, but such 
an interpretation has been shown to be untenable for the samples studied by 
the group of Grigorov [TBI EES]- Also, one of the Ioffe Institute group in more 
recent papers appears to have realised that high electric fields are not needed 
to produce highly conducting channels through polymer fihns [2~4l [231 12"5] . An- 
other member of the original Ioffe group, now in Israel, has found switching to 
high-conductivity states in channels through films of PDMS [27] . 

Baran et al. found no evidence of unusual magnetic effects in bulk atactic 
polypropylene[3S], but in this case the oxidation was probably not great enough, 
except perhaps near the surface, for conducting channels to form. I have been 
informed that a verbal report on some of the samples of Grigorov et al. was given 
in the 1990's to M. Goldes (CEO of the company with which Prof. Grigorov 
was associated at the time) by P.M. Grant, who concluded that there was no 
superconductivity. However, since no details are available, the reliability of 
this assessment cannot be judged. Dr. Grant wrote a spoof essay about the 
disovery of superconductivity at T > 541 K in polymers in 2028 in Physics 
Today [29 , but in my view put the date at about 40 years after the real event! 
Professor Grigorov appears to have stopped publication on the subject since 
2000, concentrating instead on improving samples etc., and patent work [301, 13T1 
[32] . and so, unless another group takes up studies in the field, it may be a long 
time before his work is properly assessed. 

Only one group, originally at the Institute of Synthetic Polymeric Mate- 
rials in Moscow, has consistently claimed room-temperature superconductiv- 
ity in channels through films of OAPP, PDMS and other polymers, including 
polvoctvlmethacrilate [55] . However, an early paper from the Ioffe Institute, 
and also more recent work by some scientists from the Ioffe Institute, espe- 
cially Ionov, who was a member of the original group publishing in 1990, col- 
laborating with others from Russia and Germany, have found evidence that 
superconductivity occurs in narrow channels through many polymer films at 
temperatures below the T c of superconducting electrodes, and that this is not 
due to the superconducting proximity effect. Such evidence includes conduc- 
tivity greater than 10 Q cm in channels through films between supercon- 
ducting Sn electrodes [4], Josephson-like I-V curves for superconductor-polymer- 
superconductor sandwiches in OAPP[4j, PDMS 24 , poly(phthalidylidene biphcny- 
line) [53], polvamidine [25] . and two other polymers [55]. oscillations of resis- 
tance as a function of magnetic field in films of polyimide, poly(phthalidylidcnc 
biphenylene) and PDMS below the T c of Sn[55] [231 EH], and a critical current 
for a destruction of Josephson-like I — V curves which depends on tempera- 
ture in a way inconsistent with superconductivity due to the proximity effect 
in PDMS [24] and poly(phthalidylidene biphenylene) |34 [ I37 j . There is also a 
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large number of papers by the group of A.N. Lachinov in Ufa on highly con- 
ducting channels through polymer films (see e.g.[38j[39] or a review [40]), but, 
except in recent work with Ionov and others[4Tj. without suggestions of super- 
conductivity. The work from 1998 onwards sometimes made use of spin-coated 
films, which method of production allows better control of film thickness [41] . 
Recent work of Ionov and Rentzsch[321 S3] considers the possibility of resonant 
tunnelling to explain their results. A difference between results of Ionov et al 
and those of Grigorov's group is that Ionov and coworkers report a dependence 
of conductance of channels on the work function of the substrate [U], whereas 
no such dependence is reported by Grigorov et al.. Possibly the difference is 
related to how good the contacts with the substrate are. A detailed discussion 
of different forms of / — V curves for different types of contacts to OAPP and 
PDMS films is given by Kraev et al.[45j. 

A theory for why conducting channels exist in films of elastomers (ma- 
terials with very low elastic constants) has been developed by Grigorov and 
coworkers [4H1 IT71 25J. The theory involves the formation of one-dimensional 
channels with diameters of the order of nanometers due to the formation of 
" superpolarons" , i.e. one-dimensional strings of a certain type of polaron in- 
volving electrons interacting with rotatable dipolar groups. An unusual theory 
for superconductivity in such channels involving motion of oppositely directed 
charge density waves was developed by Grigorov [351 ISO]- An earlier theory by 
me for the superconductivity involving large enhancements of pairing interac- 
tions at certain high drift velocities [51] [52] has had to be abandoned, because 
recent work on bipolarons in one dimension [53] gives no support for the greatly 
enhanced interactions at such velocities predicted by perturbation-type meth- 
ods. However, work [3] involving Bose condensation of bipolarons (without any 
special enhancement of interactions at high drift velocities) , with an initial lin- 
ear dispersion at low wave vectors may be on the right track, and the purpose 
of this paper is to improve on the model introduced in 2005, making more use 
of both experiment and theory to help to determine parameters in the theory. 

The model and constraints we use differ in four ways from that of 2005. First 
the minimum number of nanofilaments required to support superconductivity 
at room temperature is taken to be considerably higher than in the 2005 paper. 
This is because we have realised that, in the case of superconducting channels 
consisting of single-walled nanotubes between normal-metal electrodes, there is 
a minimum resistance equal to Rq/2N for singlet superconductivity, where N is 
the number of nanotubes [54] , and Rq = 12.9kf2 is the resistance quantum. For 
nanotubes there is an orbital degeneracy as well as the spin degeneracy, which 
gives rise to the factor of two in the denominator of the above expression. For 
nanofilaments with no orbital degeneracy and triplet superconductivity with all 
spins in the same direction, as is thought to be the case for OAPP, we assume 
that Rq/2N becomes 2Rq/N, where TV is the number of nanofilaments. 

The second difference from 2005 arises since we doubt if the form of disper- 
sion for bosons assumed in 2005 based on a Cooper-pair model[55] is appropriate 
in the Bose-gas regime, where the Fermi surface no longer exists even above T c . 
We assume instead that the dispersion has a Bogoliubov-type form, but with the 
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strength of the initial linear slope at T c (related to an effective bosonic chem- 
ical potential) smaller than at low T. This hypothesis is made as an attempt 
to account for the fractional increase in T c of an interacting Bose gas due to 

1 /3 

interactions, proportional to n b for low boson concentrations n& [55] [57], in 
a model of non-interacting bosons with a different type of dispersion from the 
bare one, in a similar spirit to Fermi-liquid theory. At higher boson concentra- 
tions the fractional change of T c of an interacting Bose gas due to interactions 
decreases with rth|58]. This can be taken into account in a non-interacting bo- 
son model by a mass which increases with rib, but with a dependence on rif, 
with a higher power than 1/3. Although there is no a priori reason to choose 
a Bogoliubov-type of dispersion at T c , a possible justification is that above T c 
there are fluctuating regions of superfluid which support this type of dispersion 
locally, in analogy to the Aslamozov-Larkin theory[SH] for BCS superconduc- 
tors. The steepness of the initial linear term, proportional to the square root 
of the effective bosonic chemical potential, is chosen to give approximate agree- 
ment of T c for the non-interacting bosons with a Bogoliulobov-type of disper- 
sion with numerical calculations mentioned above for an interacting Bose gas 
at low boson concentrations. Although the Bogoliubov dispersion is modified 
in a charged Bose gas (see e.g. [50]), we do not expect great changes from the 
simple Bogoliubov form at temperatures close to the condensation temperature 
because screening will cause the interactions to become short range at boson 
concentrations of interest. 

The third difference from 2005 is that we make more use of theoretical work 
(see especially [35]) to put constraints on the many parameters in the theory. We 
determine parameters by minimising a weighted sum of moduli of differences of 
two sides of equations given by constraints, with the weighting varied according 
to how accurate we think the constraints are - several of these are thought 
to be at best only of ordcr-of-magnitude accuracy. We estimate the linear 
carrier concentration for the 50 /im film from the theory of Grigorov 48 : which 
minimises the energy of a ferromagnetic nanofilament, but, in the 0.3 /im him we 
use the linear carrier concentration as an adjustable parameter because charge 
injection from the electrodes may influence the carrier concentration for films of 
thicknesses less than about 3 /im[6T] (see also Ref. [62], although the thickness 
to which this may occur does not appear to be discussed in this reference), and 
this injection may either increase or decrease the total electron concentration. 

There are twelve parameters which need to be determined in our theory. 
These are: 1) and 2) The numbers nx of nanofilaments in each transverse di- 
rection of assumed square arrays of nanofilaments for the highest-resistance 
superconducting channels for 50 and 0.3 /im films; 3) The lattice constant o,t of 
these arrays; 4) and 5) The product Cbclt of the linear boson concentration cb in 
a nanofilament with this lattice constant in 50 /im films, and the change Acbo-t 
of this product in the 0.3 /im films because of possible effects of carrier injection 
from electrodes; 6) The ratio of ar/p of ar to the radius p of the cylinder around 
a nanofilament where all dipolar groups are pointing away from the axis of the 
nanofilament; 7) and 8) The binding energy E p b and the longitudunal mass Ml 
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of the pairs which form in the polaron string due to mediation of high-frequency 
excitations; 9) The ratio p — Ml/Mt of the longitudinal and transverse masses 
Ml and My of bosons in a nanofilament array; 10) The energy of the exci- 
tation which mediates the pairing attraction between the polarons within the 
superpolaron array, which excitations we assume for most of the paper to be 
plasmons; 11) The effective bosonic chemical potential /j,b which appears in the 
assumed Bogoliubov form of the boson dispersion at T c ; and 12) The reduction 
ATi of the superconducting transition temperature for the highest-resistance 
superconducting channel for the 50 fj,m films due to phase slips. We do not 
have a method of determining the last parameter without a detailed theory of 
phase-slip resistance for wires with BEC-type superconductivity, but, assum- 
ing T c remains above room temperature for all cross sections larger than the 
minimum required to support superconductivity at that temperature, we can 
determine a lower limit on ATi « 7 K within the framework of our model. We 
do most of our calculations with ATi = 15 K, not far above the lower limit of 
about 7 K. 

After writing down the basic equations of our model in section 2, in section 3 
we find values for the first two parameters from Fig. 1 assuming good contacts. 
Then, from a combination of theory and use of experimental data, sections 2 
and 4 to 11 enable us to hnd nine equations of varying degrees of approximation 
to be satisfied by parameters 3), 4), and 6) to 11), assuming a value of 15 K for 
ATi. In section 12 we find the change of cbo-t for the thinner films to make T c 
equal to room temperature for the tit for the highest-resistance superconducting 
channel for these films. 

Some of the equations are thought to be fairly accurate while others are 
at best order-of-magnitude constraints. For our numerical work to determine 
parameters in section 13 we weight the various constraints according to our 
guesses at their accuracy. Some discussion is given in section 14. 

2. Model of Non-interacting Bosons in an Array of Nanonlaments 
with a Bogoliubov Type of Dispersion at T c 

As in Ref. [3J, we use a model of conducting channels composed of non- 
interacting bosons on a square lattice forming a square array of quasi one- 
dimensional filaments of length L and transverse dimension D measured from 
a point half a filament separation ar beyond the centers of the edge filament 
rows. The filaments are assumed to be sufficiently narrow that the bosons in a 
filament are in the lowest possible state with respect to motion perpendicular 
to the filaments at all temperatures of interest. We suppose that boson wave 
functions vanish at the edges of the array, and neglect any periodicity in the 
potential along the length of the array. Then we suppose that the boson states 
in the array can be characterised by the magnitudes of wave numbers K\ , K% 
and K3 in the two transverse and longitudinal directions, where K\ and K2 take 
on discrete values ranging from ir/D to nrin/D), and K 3 takes on values of 
kn/L, where k runs from 1 to 00. For computational purposes we approximate 
the sum to infinity by a finite sum with a suitably high upper limit. 
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The difference in our modei from that used in Ref. [5] is in the form of the Ex_ 
chosen for the bosons. Here we assume that the boson energy E B {K\, K 2 , K 3 ), 
where K\ = iir/D, K 2 — jit/D, K 3 — kn/L, has the Bogoliubov form 

E B (K U K 2 ,K 3 ) = 

[{A L {n/Lf [{P-\)+p(L/Df(P + p-2)]+^ B Y-^. (1) 



Here 
and 



A L = h 2 /2M L , (2) 

p = M L /M T , (3) 

where Ml and Mr are the boson masses in the longitudinal and transverse 
directions; fis is an effective chemical potential of the bosons at the temperature 
being considered. We say an effective bosonic chemical potential because the 
actual bosonic chemical potential at the Bose-condensation temperature must 
be zero, but, as mentioned in the Introduction, we suppose that interaction 
between bosons may be taken into account by a Bogoliubov form of dispersion 
for non-interacting bosons which give about the same type of increase in T c 
as accurate theories of the effect of interactions find. The same type of model 
with a simple quadratic dispersion was discussed in Ref. [63], but with the 
difference from the case assumed here that a discrete lattice was considered in 
the longitudinal as well as transverse directions. 

Condensation into the lowest-energy state cannot take place if the sum of 
the boson occupation factors over all states except the lowest for a chemical 
potential situated at the energy of the lowest state is greater than the total 
number of bosons, N B - Below the temperature at which this sum is equal to 
Nb , the occupation of the ground state and of states in its neighbourhood will 
start to become macroscopic, i.e. of the order of Nb, and so the equality of 
the sum with Nb determines the condensation temperature T B . For an average 
linear concentration cb of bosons per filament, this criterion can be written as 



E 



, f exp E B {Ki ,K 2 ,K 3 )/k B T B }-l 



c B Ln?p, (4) 



where the sums over K± and K 2 are sums from i, j = 1 to hit of K\ — in/ D and 
K 2 = jn/D, and the sum over K 3 is a sum from k = 1 to infinity of K 3 = kn/L, 
but the ground level i = j = k = 1 is excluded. 

As discussed in Ref. [3] , it is difficult to define unambigously a condensation 
temperature for small arrays because of the gradual way in which occupation 
numbers of the ground state change with temperature for such arrays [64j [65] . 
However, we shall be concerned here mainly with arrays containing more than 
about 10,000 bosons, and so this ambiguity should not be a problem. 



3. Constraints from Histograms of Resistances with Microcontacts 

Because of the lack of dependence of the average resistance on film thickness 
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of the lowest- resistance parts of the histograms in Fig. 1, and for other rea- 
sons discussed in the Introduction, we assume that the channels associated with 
these parts of the histograms are superconducting at room temperature. Fur- 
ther, based on previous studies[5ni [3] which found that, at least for fairly low 
cross sections of channels, the Bose condensation temperature of nanofilament 
arrays increases with the array cross section, we suppose that the supercon- 
ducting T c tends to increase as the resistance measured by the two-point probe 
method decreases. Hence we assume that the highest-resistance channels of the 
superconducting group (HRSC's) have superconducting T c 's of just about room 
temperature, i.e. we assume that 

T c ,hrsc(L) « 295^, (5) 

where L=50 /im or 0.3 /im for the highest-resistance channels of the low- 
resistance groups in Figs 1(a) and (b) respectively. 

For an array of single-walled superconducting carbon nanotubes with good 
contacts from normal metals, there is a minimum resistance [54] 

R m = R Q /2N, (6) 

where 

R Q = 12.9kfi (7) 

is a resistance quantum. Single-walled carbon nanotubes have two orbitally- 
degenerate and spin-degenerate subbands[54] [66] . We assume that this mini- 
mum resistance becomes 2Rq /N for triplet superconductivity in nanofilaments 
with all the spins of a nanonlament in the same direction, as is thought to 
be the case for OAPP (46] |47l [48]. We have no reason to assume orbitally- 
degenerate subbands, and, if all spins are aligned, this also halves the number 
of current paths. Hence, assuming good contacts, we suppose that for the 
highest-resistance superconducting channels for the two film thicknesses being 
considered satisfies 

n T ,HRSc(L) 2 « 2Rq/R, (8) 

where R is the resistance for the relevant channel. Using data from Fig. 1 we 
deduce 

riT,HRSc(0.3fim) w 17, ht,hrsc (50/xm) sa 32. (9) 

For samples with small cross sections, phase slips cause resistance [671 [68] . 
Defining, arbitrarily, T c of the system as that temperature at which phase slips 
cause a resistance of half the resistance in the normal state, then phase slips 
reduce an initial superconducting transition temperature T c0 to T c , with 

T c = T c0 - AT (10) 

where, for BCS-type superconductors [67], 

AT oc (cross — section) -2 ^ 3 oc n T 4 ^ 3 . (11) 
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Although there are papers on phase slips for Bose condensates [521 [ZD] > they do 
not appear to give sufficient information to determine the dependence of AT 
on fix or on coherence lengths for such systems. Thus we use Eq. (11) to find 
the dependence of AT on rip for a given boson concentration in our numerical 
work. 

4. Relation between Linear Boson Concentration and Other Param- 
eters based on Grigorov's 1991 theory 

In an elastomer (a material with very low elastic moduli) containing rotatable 
groups of atoms with electric dipole moments, an excess electron causes all 
dipoles to point away from the electron within a radius Rq given bv|46j 

R = (eD m /TTvGe h )i, (12) 

forming an unusual type of polaron. Here D m is the dipole moment of a group, 
v is its volume, G is a shear modulus, and tu is a high-frequency dielectric 
constant. For many electrons in such a material it is energetically favourable 
for the polarons to form strings called " superpolarons" [151 271 SB] • For super- 
polarons with a linear concentration of polarons cp , the dipoles point outward 
in a cylinder of radius p = 2RqCf^8\. In our theory we suppose that there is 
pairing of polarons within a string caused by attraction mediated by plasmons 
or high-energy phonons, to give a linear concentration of bosons cb = 0.5cf, 
and so, putting results in terms of boson concentrations we have 

p = mlc B . (13) 

When exchange forces are included, Grigorov[48] finds that, for low carrier con- 
centrations, the minimum energy of a polaron string occurs when the polarons 
are ferromagnetically aligned with a linear concentration such that 

k F R « 1.6, (14) 

where k F , the Fermi wave number for spin-aligned electrons, is 

kp = ircp = 2itcb- (15) 

Ferromagnetism at low concentrations in a quasi one-dimensional electron gas 
has also been discussed by Shelykh et al.|71]. 

We consider a square array of polaron strings (to be referred to as nanofil- 
aments from now on) with lattice constant ap. From Eqs. (13-15), we find 
that 

c B a T = 4 x (l.6/27r) 2 (a T /p) = 0.259(a T /p). (16) 



5. Relation between Separation of Nanofilaments and Radius of Fully 
Aligned Dipole Region from Experimentally Based Estimates of Av- 
erage Carrier Concentration in Conducting Channels 
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The average concentration of bosons within channels is not known accurately. 
However, we may obtain an order-of magnitude estimate of this as follows. From 
magnetisation results, the maximum concentration of aligned spins in any sam- 
ple reported on up till now appears to be that for sample 1 of Ref . JSHj , where 
the saturation magnetisation after any diamagnetic contribution has been re- 
moved by field cycling is 0.058 emu cm -3 (see dotted curve in Fig. 1 of Ref. 
[20]). Hence, assuming this moment is due to aligned electron spins, the average 
concentraion «f of fcrmions in this sample is given by 

n F « 6.3 x 10 18 cm" 3 . (17) 

The highest concentration of conducting channels reported in any sample is 
given in Ref. |72] , where the average separation of the channels was determined 
to be 7.5 /im. Hence with a typical channel diameter of about [5] 1 /xm, the ratio 
V s /V c of the volume of the sample to the total volume occupied by conducting 
channels is given by 

Vs/Vc ~ 7.5 2 x (4/tt) = 72. (18) 

If we suppose that the fraction of material occupied by channels in sample 1 of 
Ref. [20] is about the same as in the sample of Ref. [72], then, from Eqs. (17) 
and (18), we deduce that the concentration of bosons ub in the channels (half 
the fermion concentration) satisfies 

n B ~ 2.2 x 10 2O cm" 3 . (19) 

Using 

riB = (c B a T )a^ 3 (20) 
and Eqs. (16) and (19), we deduce that 

a T « 1.04(a T /p) 1/3 nm. (21) 

Although we used ~ in Eq. (19), we have changed to « in Eq. (21) because of 
the one-third power. 

6. Estimate of Transverse Bipolaron Mass based on Grigorov's 1991 
Theory 

In this section we make use of the theory of Ref. [48] to estimate the transverse 
mass for motion of particles across our array before inclusion of pairing due to 
mediation of plasmons or high-frequency phonons, and use these calculations 
to estimate the ratio Ml /Mt of longitudinal and transverse bipolaron masses 
in terms of other parameters. Assuming a tight-binding model for transverse 
motion across the array, the transverse single-particle masses can be found in 
terms of energy overlap integrals t. 

We suppose that the diameter D of the neutral cylinder enclosing a polaron 
string in the theory of Ref. [35] can be identified with or for our array of strings. 
The potential energy appearing in an expression for t has to be measured from 
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a zero equal to the lowest energy state of an isolated string. This energy, which 
we call Eo, is given by 

/•OO 

E = 2tt {x 2 (r)V(r) - (h 2 /2m e ) X d 2 X /dr 2 }dr, (22) 
Jo 

where the (real) radial wave function xi r ) is given by 

X (r) = (l/2^) 1 / 2 (2/3)2 1 /2 /3 (i + /3r)cxp(-/3r), (23) 

with 

$ = 5/R = 5 x 2ttc b /1.6 = 19.6c B (24) 

in our notation. We have changed slightly from Grigorov's notation in that we 
have included a factor (1/2tt) 1 ^ 2 in \- The potential energy of an electron V{r) 
in our notation satisfies 

V(r) = -[e 2 (l-k)/(e h a T )} 
{2c B a T [2\n{a T /2p) + {2p/a T ) 2 - 1] + a T {p - r)/i? 2 } (r < p) 
= -[2e 2 (l - k)c B a T /(e h a T )][2Ha T /2r) + (2r/a T ) 2 - 1] (p < r < a T /2) 

= (a r /2 < r). (25) 

Here 

fc = e;,./e s , (26) 

where e s is the static dielectric constant, and eh, is the high-frequency dielectric 
constant. For materials in which polaron strings occur, k is small compared 
with unity. 

For wave functions as given by Eq. (23), we find that the energy overlap 
integral t is given by 

/■OO />27T 

t= / 2Trr[V(r) - E )]x(r)x(u)drdd, (27) 
Jo Jo 

with 

u = [a% + r 2 - 2a T rcos6] 1/2 . (28) 

The bare transverse mass m FTb of fermions for our square lattice in a tight- 
binding model is related to the energy overlap integral t by 

m F Tb = h 2 /2ta 2 T . (29) 

Here we use bare in the sense of not including any effects of interactions with 
high-frequency modes which give rise to the pairing in our model. 

For a Holstein model with anisotropy[73], the ratio of the transverse to lon- 
gitudinal polaron masses is not much different from the same ratio for the bare 
masses, at least for interactions with high-frequency phonons. Thus, combining 
this result with the assumption that the ratio of the pair transverse mass Mr 
to the bare transverse fermion mass m ft is approximately equal to the ratio of 
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the pair longitudinal mass Ml to the bare fermion longitudinal mass, which we 
take as equal to the free-electron mass m e , we find 

M T « 2m FT b(M L /2m e ). (30) 



7. Estimate of Energy of Mediating Bosons assuming these are Plas- 
mons 

We assume for most of this paper that the bosons mediating the attraction be- 
tween fermions in the nanofilaments are plasmons, although we make some 
remarks in the discussion section about the alternative that they are high- 
frequency phonons. 

As in Ref . [5T] , we estimate the plasmon energy as that for a single nanofil- 
ament at wave numbers q such that qd ~ 1, where d is the diameter of the 
nanofilament. In this case the plasmon energy is 

hcu = h(mire 2 c B d- 2 /M L e h ) 1/2 , (31) 

where cb is the linear boson concentration and th is the high-frequency dielectric 
constant. The difference from the usual factor of Air in the bracket comes about 
since the pair charge is 2e. For the diameter d in Eq. (31), we first estimate the 
mean value of the radius from the axis from Grigorov's variational method [15] 

as 

/>oo y* oo 

r m = / r 2 (l + Pr) 2 e- 2tir / / r(l + ^rfer 2 ^ = (14/9)/r\ (32) 
Jo Jo 

where, from Eq. (24), /3 — 19. 6 {c b <Xt) /clt- For a uniform distribution of charge 
up to a cylinder radius of r c , the mean value r m of r is given by 

r m = (2/3)r c . (33) 

Thus to convert the mean value of r given above to a total cross section, we 
multiply the cross section by (9/4), or the linear dimensions by (3/2). Hence 
the diameter of an equivalent square filament is taken as 

d = (3/2)(14tt 1 / 2 /9/3) = 0.2lla T /{c B a T ). (34) 

Hence, from Eq. (31), 

fiuj w h(167re 2 CBaT/e h M L a^) 1/2 (4.7AcBa T ). (35) 

There will be additional terms proportional to the square of the wave number 
and inversely proportional to the mass at high q, but, for parameters of interest, 
these are not expected to change the relevant average plasmon energy by large 
amounts, and will be partially cancelled by the reduction in plasmon energy at 
small q. 
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8. Approximate Relation between the Bipolaronic Binding Energy 
and the Bipolaron Mass 

We obtain another approximate relation between parameters by making use 
of results for bipolaron binding energies and masses for weak to intermediate 
coupling in a one-dimensional system 53 . Although these results were obtained 
with symmetric spatial bipolaron wave functions, implying singlet bipolarons, 
results for singlet and triplet bipolarons in a Holstein-Hubbard model do not 
look drastically different from each other for weak to intermediate coupling (see 
Figs. 1 and 2 of Ref. [7J or Fig. 15 of Ref. [75]). From section 4 of Ref. [5J, 
for V 2 — 0.125, 0.25 and 0.375 and a = 4 in the notation of that paper, a = 4 
corresponding to high frequencies of the bosons mediating the electron-electron 
attraction, we can see that the ratio (Eb/fou))/[(Mb p /2m,b) — 1] lies between 1.63 
and 1.32, and, with extrapolation to higher V 2 a up to 0.5, this value could go 
down to about 1.25. For values of V 2 between 0.125 and 0.375, the values of 
the bipolaron binding energies E p b lie between 0.08?ki; and 0.387k*;, while the 
bipolaron masses Ml (for a bare mass mj) lie between 2.13 and 2.5 nib for the 
same parameter range. Based on these results we shall suppose that 

M L /m b = 2.0+ (lA±0.2)(E pb /hu)), (36) 

and in our case we suppose mj, = m e , the free-electron mass. 

V and a in our 2007 model, which uses local electron-boson interactions but 
a quadratic dispersion for bare electrons, correspond approximately to what are 
usually called g and {hui/i) 1 / 2 in the Holstein model. Thus Q.bV 2 a 2 corresponds 
to g 2 huj/2t, usually called A in the 1-D Holstein model. Hence V 2 a = 0.5 and 
a = 4 imply A ps 2. Fig. 1 of a paper by Hohenadler and von der Linden [75] 
shows that the kinetic energy (related to reciprocal masses) of Holstein bipo- 
larons is reduced by a factor of more than five at A=2. However, we expect 
reductions of masses in our model compared with those on the Holstein model 
to become larger as the coupling increases, and so the comparatively small 
masses for the values of V 2 and a with which we have been concerned may be 
realistic for our model. 

9. Probable Lower and Upper Limits for the Bipolaron Binding En- 
ergy 

Although there are criteria for being on the BEC side of the BEC to BCS 
transition in terms of the product of the scattering length and the Fermi wave 
vector [771 1751 [75] or in terms of the ratio of E p b to Ep (deducible from Fig. 1 of 
Ref. [80]), we are not sure how these are modified for systems which are both 
polarised and in which there is also a cut off in wave vector in two directions, 
with the Fermi wave vector having a value which is approximately equal to its 
1-D value. (How good this approximation is will be discussed after numerical 
work determined approximate values of parameters). 

The simplest criterion for being on the BEC side of the transition is that 
the diameter of the pairs is smaller than their mean separation. In the longi- 
tudinal direction, the diameter of the pairs is approximately {Sfi 2 /MlE^) 1 ^ , 
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assuming Ml k, (twice the longitudinal fermion mass), and the pair separation 
is approximately (cbOt) -1 ot ■ Hence a lower limit for the pair binding energy 
is 

E pb = 16{h 2 /2M L a 2 T )(c B a T ) 2 . (37) 

Another lower limit on the pair binding energy requires that E p b is appre- 
ciably larger than 3.5 k B x 295 K « 0.08 eV (BCS value of energy gap for 
T c = 295 A'). We use the first of these limits, but after numerical work we also 
check that the value found is also above the second lower limit. 

For an upper limit, based on our discussion in the previous section, and the 
assumption that we do not wish to venture far out of the intermediate-coupling 
regime, we suppose that the pair binding energy E pb satisfies E p b < 0.6hu>. 
Hence, using Eq. (37), we put 

E pb ~ 0.5[16?i 2 /{2M L a 2 T ){c B a T ) 2 + OMlo]. (38) 



10. Order of Magnitude Constraint from Critical Current Density 

Estimates of the critical current density made from pulsed-current measure- 
ments vary from[5T| 10 8 Acm -2 to[S] 5 x 10 9 Acm~ 2 , depending on whether the 
size of the pressure contacts or the estimated cross section of a single channel is 
used to convert a critical current of (63 ±17) A, for gradually increasing magni- 
tude of the currents in the pulses, to a current density. We take the geometric 
mean of these two types of estimates, i.e. 

j c ~ 7 x 10 8 AcuT 2 . (39) 

In numerical studies of a fermion superfluid, Heiselberg[82 shows that, for 
the Bose-gas regime, the sound velocity at low temperatures (equal to the critical 
drift velocity Vd for a Bose gas) lies between about 0.08 and 0.25 times the Fermi 
velocity, vp. Thus we assume that 

v d = (0.165 ±0.085)i; F . (40) 

Since the critical temperature in the pulsed currents used is estimated in Ref. 
[5] to be greater than 700 K, we assume that room temperature will count as a 
fairly low temperature for the purposes of estimating the critical current. Hence, 
using Eqs. (39) and (40), we find for our model that 

(0.33 ± 0.l7)e(c B a T )a T 3 v F - 7 x 10 8 Acm~ 2 = 2.1 x 10 18 esu. (41) 

We approximate the Fermi velocity by that for a one-dimensional ferromagnetic 
system as in the previous section. Then we find 

v F w h(2Tr)c B /m FL w h(4ir)(c B a T )/a T M L . (42) 

Eqs. (41) and (42) give another order-of-magnitude relation between parameters 
of our model. 
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11. Estimate of Effective Pair Chemical Potential 

We estimate an effective pair chemical potential to put into the Bogoliubov- 
type of dispersion for non-interacting bosons at T c as follows. First we look at 
accurate results for the condensation temperature for interacting bosons [561157] 
at low carrier concentrations, where it is found that the change AT C from the 
non-interacting boson value T c0 is given by 

AT C /T c0 « l.3a B nJ\ (43) 

where as is the boson scattering length. Equation (43) is only valid for low bo- 
son concentrations, but decreases of T c due to interactions for higher concentrations^ 
can be taken into account in a non-interacting boson model by assuming an ap- 
propriate increase in mass with n b with a dependence on density with a higher 
power of riB than 1/3. 

The boson scattering length in Eq. (43) satisfies 

a B « 0.6a F , (44) 

where op is the fermion scattering length|83 , which in turn is related to the 
pair binding energy E p i, and fermion mass rap by [84 

h 2 /(m F a 2 F ) = E pb . (45) 

Assuming the fermion mass is approximately equal to half the boson mass 
Mb, we deduce from Eqs. (44) and (45) that 

a B « 0.6(4?i 2 ' l2M B E ph ) x l 2 . (46) 

For an anisotropic system we suppose that as to use in (43) is as in (46), but 
with 

1/M'J 2 -+ (1/3)(1 + 2p 1 /2 )/M V2 (47) 
in our notation. Also, in our system, 

n 1 ^ = (c B aT) 1/3 a^\ (48) 

Hence, from Eqs. (43) to (48) we find 

AT C /T c0 « 0.52(^ 2 /2Ml^4) 1/2 (1 + 2p 1 / 2 )(c s a T ) 1/3 (49) 

For non-interacting bosons with a Bogoliubov type of dispersion at T c , we 
find numerically the [Ib which gives the same fractional change of T c from that 
for a pure quadratic dispersion as for the interacting system as determined 
above. Hence we find \ib as a function of our other parameters. 

12. Change in Carrier Concentration in 0.3 /im Films required to 
make the T c of the Highest-Resistance Superconducting Channel for 
these Films equal to Room Temperature 
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Now, because of possible injection of carriers from the electrodes (see the Intro- 
duction), we assume that 

c B a T -)• c B a T + A(c B a T ) = c B2 a T (say) (50) 

for 0.3 fim films. Then p in our equations in section 2 has to be replaced by 
P2 = Ml /Mt2 , where the transverse mass Mt2 for the thinner films is calculated 
as in section 6, but with c b clt replaced by c B 2<it- 

For a fixed boson concentration, the decrease AT of T c is taken to be given by 
AT oc 7i T 4 ^ 3 [see Eq. (11)], but we also suppose that AT is proportional to the 
inverse of the carrier concentration in order to take into account that phase slips 
are not so easily formed if the carrier concentration rises. Hence we suppose that 
the reduction in T c , AT2, for the highest resistance superconducting channel in 
the 0.3 /im films, with tit — 17, is related to the corresponding quantity, ATi 
for 50 /im films, with ut — 32, by 

AT 2 » 2.3(c B a T /c B 2a T )AT 1 , (51) 

where cg2 is the linear carrier concentration in the 0.3 /im films. 
After other parameters are determined, c B 2Clt is varied so that 

T B - AT 2 = 295A" (52) 

for tit = 17. The change in c B ar also gives rise to a change in p — Ml/Mt 
because it changes the potential involved in calculating the transverse mass. 

13. Numerical Results 

For our numerical work we first choose a value for the parameter A^ . T B in our 
model as a function of ut reaches a shallow maximum as a function of tit , and 
then declines slightly as ut decreases further. (A qualitative discussion of why 
this slight decline occurs has been given for bosons with a quadratic dispersion 
in Ref. [63]). However, the results shown in Fig. 1 indicate that once T c reaches 
room temperature it remains above that temperature as Ut increases further, 
at least up to very high tit- Hence, when other parameters are approximately 
determined, we may obtain a lower limit for ATi from this condition. For our 
final values of parameters we find that ATi > IK, and for our calculations we 
take a value 

ATx = 15K. (53) 

That ATi is probably not much larger than this may be argued by comparison 
with conclusions of Zhao for carbon nanotubes[5J [9]. He gives an example (see 
pp 46-48 of Ref. [5]) where phase slips give rise to a resistance equal to about 
0.4 x (the normal-state resistance at a temperature of the order of 0.8 x the 
transition temperature without phase slips) for a number of (spin-degenerate) 
channels equal to 54, corresponding to Nj. ~ 108 for our type of system, where 
we suppose there is no spin degeneracy. Hence, if we assume that Eq. (11) is 
true, for our case we might suppose that AT\/T B for 32 2 = 1024 nanofilaments 
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is slightly larger than [(1024/108)" 2 / 3 = 0.223] x 0.2 w 0.045, implying AXi, 
given by A7\ = 0.045 x (295K + ATi) is slightly larger than 14 K. Of course 
this is only an order-of- magnitude estimate because Eq. (11) and the analysis 
of Zhao are based on a theory appropriate for materials on the BCS side of the 
BCS-BEC transition, whereas we are assuming we are on the opposite side. 

Having determined tit for the highest-resistance superconducting channels 
in section 3, and chosen the value of ATi given by Eq. (53), we next make use 
of the requirement that T c s=y 295K for the highest-resistance superconducting 
channels for 50 /jm films and results of sections 4 to 11 to determine another 
eight parameters of our theory approximately. These are ar, cbclt for 50 /zm 
films, a^/p, Eb, Ml, p = Ml/Mt, and [1b- This gives nine approximate 
equations for eight parameters, and so one might think that the parameters are 
overdetermined. However, because many of the relations between parameters 
are very approximate, we do not have to worry about satisfying them exactly, 
and it turns out that we can satisfy them to within our estimated accuracy for 
the equations. 

Our procedure is as follows. We regard Eq. (16) of section 4 as exact, 
and then minimise a weighted sum of moduli of differences between sides of 
dimcnsionlcss quantities based on other equations, with weightings based on 
our estimates of accuracy of our equations. In order not to cause computing 
time to be too long, we approximate the sums required to produce Tb given in 
section 2 by grouping terms in the sums over k in section 2 as equal to their 
mean values between two limits, the differences between these limits increasing 
as k increases. For our detailed calculations, we regard the numerical equation 
giving T c = 295 K for the highest-resistance superconducting channel for 50 /zm 
films as accurate to less than 1%, but supplemented by errors of up to about 1% 
because of our approximate methods of evaluating the sums over k in equations 
of section 2. We weight the differences in moduli of calculated and fitted values 
of quantities x divided by the fitted values, i.e. 

\{x c -x)/x\ (54) 

by weightings according to our guesses of the accuracy of the equations. Tak- 
ing the weighting for |[(T CC - 295X)/295if]|, where T cc is the calculated T c for 
the highest-resistance superconducting channel for 50 /im films, as 1, we choose 
the following much smaller weightings for similar quantities associated with ay 
given by Eq. (21) (0.06), j c given by Eqs. (41) and (42) (0.015), p = M L /M T 
given by the methods described in section 6 (0.03), Uui based on Eq. (35) (0.1), 
E pb based on Eq. (38) (0.03), M L based on Eq. (36) (0.2), and based on 
methods described in section 11 (0.1). There is some guessing as to the appro- 
priate weightings. The weighting for the transverse bipolaron mass is taken as 
fairly low because its calculation depends both on details of Grigorov's varia- 
tional method, and assumptions about how to convert the masses deduced from 
his theory to masses for pairs. The very small weighting for j c is due to the un- 
certainty in the experimental values of critical currents, the larger uncertainty 
arising from conversion to a critical current density, the uncertainty to what 
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extent 295 K represents a "low T" for calculation of critical currents, and the 
uncertainty about how much other parameters change when large currents are 
passed. The weighting for E pb is small because for typical values of parameters 
the ratio of the upper and lower limits is quite large. For computational pur- 
poses, for calculations of p we approximated results for logarithms of expressions 
as polynomials in relevant parameters, and for /is we approximated by a differ- 
ent simple expression in the relevant parameters, and used the polynomial forms 
or other simple expressions in our minimisation programs. Since clt/p > 2, but 
probably not much greater, we only did calculations for cbclt (=0.259aT//?) 
between 0.5 and 0.7. For the eight-parameter minimisation, we made use of a 
program AMOEBA from a book[5S], based on a method suggested by Nelder 
and Mead [86]. For our computing we used dimensionless lengths in units of 
nm, dimensionless masses in units of m e , and dimensionless energies in units of 
(h 2 /2M L )(n/a T ) 2 . 

From the best fit to the weighted sum of moduli as given by Eq. (54), we 
find that 

T c — 295.03 K, a T /p = 2.061, (55) 

implying with Eq. (16) that cbo-t — 0.534, and fitted and calculated values x 
and x c of other quantities, and of [{x c — x)/x] found for (i) ax from Eq. (21), 
(ii) p from section 6, (hi) the plasmon energy from Eq. (35), (iv) Ml from Eq. 
(36), E p b from Eq. (38), (v) j c from Eqs. (41) and (42), (vi) \xb from section 
11 are: 

ot = 1.282nm, (ot c — a>T)/a>T = 0.033 

p= 0.315, (p c -p)/p = 0.002 

riw = 1.31eV {Tilo c - Tiuj)/huj = -1.3 x 10~ 6 
M L = 2.376m e , {M Lc - M L )/M L = 0.029 
E pb = 0.415eV, (E pbc - E pb )/E pb = 7 x lO" 5 
j c = 3.43 x 10 8 Acm- 2 , {j cc - j c )/j c = -0.510 

li B = 0.00127eV, {hbc - Vb)/pb = 1.1 X lO" 5 . (56) 

The only quantity for which the calculated values and the values obtained by 
our fitting process is significantly different is j c . Some reasons why we expected 
that the j c estimate would be inaccurate are stated earlier in the section. For 
the values of parameters found, the upper limit for E pb in section 9 is 0.79 eV 
and the lower limit is 0.045 eV. It turns out that our lower limit is smaller than 
the alternative lower limit of 0.08 eV given in section 9, but if we had used 
this alternative in our computation we do not expect that there would be much 
change in parameter values. 

Note that the small differences in most cases between calculated and fit- 
ted parameters shown above do not imply that the parameters are known to 
anything like that accuracy, as errors cannot be expected to be less than our 
estimated errors in the calculated relations used in our fitting process. 

After the other parameters are obtained, we treat CB20-T for the channels 
through 0.3 /im films as adjustable, where cbi is the linear concentration of 
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bosons in these channels. The calculation of transverse mass in this calculation 
has to be modified for the different carrier concentration. We find that Cbi<it = 
0.636 (compared with cbo-t = 0.534 for 50 fim films), and the the ratio P2 of 
longitudiinal to transverse mass in the channels through the thinner films is 
given by P2 — 0.258, smaller than the value p = 0.315 for the thicker films. 
Note that, if other parameters were to remain fixed, then Tb would increase 
with boson concentration. However, p = Ml/Mt decreases with increase of cb 
because of the increase in depth of the potential of Eq. (25) with increasing eg. 
Thus there is a subtle competition between two effects in the parameter range 
in which we are interested. 

For the values of parameters found, T c = Tb — 15 K x(tit/32)~ 4 / 3 reaches 
a maximum of about 303 K for 50 /im films for ut ~ 140, and T c = Tb — 15 
Kx(n T /32)- 4 / 3 (c B a T /cB2aT) for 0.3 /im films reaches a maximum of about 
310 K for tit ~ 90, assuming channels are single domains with all spins aligned. 
For relatively large channels it is likely that the channels will split up into 
domains with opposite directions of spins, but we do not know enough about 
the transverse exchange constants or magnetic anisotropy to be able to predict 
the maximum diameter of single-domain channels. In order to explain higher 
T c 's than the two values noted above, observed when relatively high currents are 
applied 3], we can invoke decreases in transverse masses due to reduced lattice 
constants produced by current-current interactions between nanofilaments as 
discussed in some detail in Ref. [3]. 

14. Discussion 

In this section we discuss many of the assumptions that come into our model, 
and to what extent they can be justified. We also discuss some other matters. 

Assumption of an array of nanofilaments 

There is at present no experimental evidence that the conducting channels are 
composed of an array of nanofilaments, and so for this assumption we have 
to rely on the plausibility of the theory of Refs. (46] [47] [48] for strings of an 
unusual type of polaron occurring in elastomers with rotatable dipolar groups. 
In this connection we remark that there are also theories for strings of more 
conventional polarons (Ref. [87] and references therein, also Ref. [88] ) , although 
to form these one has to include interactions other than the usual Frohlich 
interactions with polar phonons and Coulomb repulsions, and there are limits 
to the lengths of the strings [88 . 

Assumption of a square lattice 

We cannot rule out a triangular lattice. Our use of a square lattice is based on 
the simplicity of calculations, and the probability that results for a triangular 
lattice will not be much different, at least for weak or intermediate coupling of 
electrons and bosons mediating the attraction. For strong coupling a triangu- 
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lar lattice is likely to be favoured because of the smaller transverse masses of 
bipolarons on such a lattice when coupling becomes strong 89 . 



Assumption that our system is on the BEC side of the BEC-BCS transition 
We estimate the pair diameter d p in the longitudinal direction from the pair 
binding energy and reduced mass of the particles making up the pair, and com- 
pare this with the separation between pairs s ss (cbo-t) _1 ot- We use the 
approximation sign rather than the equality for s because there will be a slight 
deviation from this value towards the (lower) average separation [cbclt) -1 ^ 3 ^ 
which would occur if the system were three-dimensional with no cut off in wave 
vector in the transverse direction. We find d p = A(Ti 2 /2MlE p i j ) 1 ^ 2 = 0.79 nm, 
s » 2.40 nm, and so the diameter of the pairs is considerably smaller than their 
mean separation. This shows that we are well on the BEC side of the transition. 



Assumption of good contacts 

If the contacts were not good, then jit for the highest-resistance superconduct- 
ing channels would have to be larger for both thicknesses of films, and then it 
would be more difficult to justify a ATi as large as we need to permit T c to 
increase monotonically with ut- Also, from Fig. 1, the ratio of resistances of 
the lowest-resistance mediumly conducting channels imply, taking into account 
the different channel lengths, that the ratios of cross sections for the lowest- 
resistance normal channels for the two film thicknesses needs to be about four, 
the same as inferred for the highest-resistance superconducting channels assum- 
ing good contacts. Hence, if the contacts are not good, they would have to have 
about the same departure from goodness for both thicknesses of films, which 
would be another assumption. For discussion of normal channel resistances, 
although with smaller values of than given by the present model, see Ref. 



Assumption of a Bogoliubov-type of dispersion for bosons at T c 
A possible justification for this assumption is discussed in the Introduction. We 
tried to fit the data assuming a simple quadratic dispersion, but were unable 
to do this with plausible values of parameters. However, we could fit the data 
with a fixed value of /ib [in units of E u = (h 2 /2M l)(it / clt) 2 ] over three orders 
of magnitude smaller than that found by our data fitting. For hb/E u = 10~ 5 
[cf. hb/E u = 0.0132 for the results given in section 13], ay is reduced to 1.037 
nm, p is increased to 0.423, M L = 2.444m e , and a T /p = 2.000. For ^b/E u 
significantly smaller, Tg for ~ 32 is too small for plausible values of aj and 
boson concentrations, e.g. for \xb = 0, L=50 /im, a? — 1 nm, p — 0.5, ca^ — 
0.54, Ml = 2.55m e , Tg for =32 is only 25K, compared with nearly 500K 
for large tit- The only possibility to fit data with /is = would be to assume 
poor contacts, so that the lowest n-r to give T c — 295K is considerably larger. 
We tried fitting assuming that nr for the highest-resistance superconducting 
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channel for the 50 /im films was nj- = 100 (setting AT = in this case because 
of the large Ht), but did not get nearly so good a fit as for our Bogoliubov-type 
dispersion and tit = 32. For kit = 100, we found that the calculated value of ax 
was 51% higher than the fitted value, the calculated value of p was 72% lower 
than the fitted value, and the weighted sum of moduli of differences between 
calculated and fitted values was about seven times higher than for our previously 
reported results for ut = 32 and fiB as a fitting parameter. The problem with 
trying to fit with low values of tit for /is = is that, for quadratic dispersion, 
the rise of Tb with tit is rather slow, and so, if we use Grigorov's expressions 
for linear carrier concentrations cb, it is not possible to get Tg for small tit 
up to room temperature unless ot is small, and then it becomes impossible to 
satify Eq. (21), since we have the restriction that ot/p > 2. Another problem 
is the strong length dependence of Tb for quadratic dispersion (proportional 
to 1/L for small tit), and so one would need to have a considerable change in 
carrier concentration in the 0.3 /im films to find a room-temperature Tb for 
the highest-resistance superconducting channel with parameters determined by 
fitting Tb for the highest-resistance superconducting channel for 50 /im films. 
Since computation time for our program goes up as the square of Ht, it is not 
practicable to see whether things improve for pure quadratic dispersion if tit 
is taken as much larger than 100, without making further approximations for 
Tb- For the rather drastic approximation made by Eq. (21) of Ref. [53], but 
increased in magnitude by 15%, since, for the cases we studied there gave errors 
of the order of 15% in Tb, we found that the best fit was obtained when tit for 
the highest-resistance superconducting channel for the 50 /im film was about 
240, with values of other parameters as follows: 

p = 0.596, clt/p = 2.003, a T = 0.942 nm, M L /m e = 2.444, huj = 1.96 eV, 
Eb = 0.621 eV, j c = 1.08 x 10 9 A cm -2 . However, neither the goodness of the 
fit nor values of all parameters vary monotonically with tit, which may indicate 
that we do not always reach the absolute minimum of our weighted sums with 
our program. For these values of parameters, tit^t for the highest resistance 
superconducting channel for 50 /im films is approximately equal to 230 nm. 
Since there is a spread of resistance for the low-resistance channels for these 
films shown in Fig. 1 of about 100, this would imply that the widest channels 
shown in Fig. 1 would be about 2.3 /im in diameter, slightly larger than the 
largest possible value mentioned in Ref. [5]. 

We may look at results for pure quadratic dispersion in more detail some- 
time, if we can find a better approximation to Tb which does not require too 
much computing time to investigate. Although we think that our arguments 
about fluctuations giving parts of the samples at any one time with higher Tg's 
than average, and hence a Bogoliubov-type of dispersion for these parts is plau- 
sible, it is less clear whether it is appropriate to consider some average type of 
Bogoliubov dispersion over the whole sample with a lower effective value of [Ib, 
as we have done for most of this paper. 

Remarks on Grigorov 's theory for ferromagnetic nanofilaments 
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It might be objected that a paper by Lieb and MattisJSD] could rule out ferro- 
magnetism in quasi 1-D filaments. However, on studying this paper, it appeared 
to the present author that Lieb and Mattis's result depended on the system be- 
ing really one dimensional, not just quasi one dimensional in the sense that all 
the carriers are in the lowest state with respect to transverse motion, and thus 
we think their result is not relevant to any real system. Grigorov[48j also notes 
that the result of Lieb and Mattis is based on a model with short-range forces. 
Further, Shclykh et al.|71) also find ferromagnetism in a quasi one-dimensional 
electron gas, and Meyer and Matveev|91j mention, in connection with quantum 
wires, that "deviation from true one dimensionality may, in principle, give rise 
to a spin-polarised ground state of the interacting electron system" . 

There are more general objections to any form of long range order in quasi 
one-dimensional systems, but these will not be relevant for arrays of nanofila- 
ments. 

Stability of electron strings 

Grigorov[25] finds that the ferromagnetic electron string is only stable if i?o, 
the radius from an excess electron up to which all dipoles point away from the 
electron, satisfies Rq > 0.6 nm. For our parameters and equations in section 4 we 
find i? ~ 0.61 nm. Putting this value into Grigorov's Eq. (9), with the second 
term reduced by a factor of 0.8 since we are assuming k = e/j/e s = 0.8 instead of 
unity as in his calculations, and using — 2.2, we find E min — (0.96 — 0.61) = 
0.35 eV. However, in our model there is a reduction in energy due to polaron 
and bipolaron binding energies caused by interaction with high-energy bosons. 
From our fitting we find that the bipolaron binding energy E p b ~ 0.42 eV, 
with E p b/hoj=0.32. From Ref. [53] . section 4, for a = 4, corresponding to 
fruj/t = 16, this value of E p b/hui probably corresponds to V 2 « 0.29. Then for 
weak coupling and a — 4, from equation (46) of Ref. [52], the single-polaron 
binding energy E s b = 0.28huj = 0.37 eV. Hence we get a total energy lowering 
per electron of E s b + 0.5E p b — 0.58 eV, quite enough to stabilise the system. 

Note that the value of i?o = 0.61 nm noted above that we find is considerably 
smaller than the estimate that i?o is greater than about 2 nm mentioned in Ref. 
|48j . and so, if our parameters are approximately correct, then we think Grigorov 
must have underestimated the shear modulus G in Eq. (12). Other quantities 
in this equation are probably known far better than G. 

Alternative assumption for Hlu 

There are high-frequency phonons in polypropylene with energies Tiuj — 0.37 — 
0.39 eV 92 . We have also looked at our model on the assumption that phonons 
of energy 0.38 eV mediate the attraction. The best fit when minimising the 
weighted sum of the same quantities as before except for the calculated and 
fitted frequency of bosons, gives ax/p = 2.055, and the following values of other 
parameters and departures from their calculated values: 

ot = 1.322nm, (ot c — ar)/aT = — 8 x 10~ 7 
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p = 0.234, (p c -p)/p = 9x 1Q- J 

M L = 2.527m e (M Lc - M L )/M L = -7 x 1(T 4 

E pb = 0.143eV, {E pbc - E pb )/E pb = -0.063 

j c = 2.84 x 10 s Acm- 2 ( Jcc - j c )/j c = -0.595 

/ib — 0.00263eV, {^ Bc - Vb)/vb = 1.1 x 10~ 7 . (57) 

Differences of calculated and fitted parameters are on average not much 
different from values obtained assuming plasmons mediate the interaction. For 
these parameters we find d p = 1.29 nm and s rs 2.48 nm, and so we are still on 
the BEC side of the transition, although not quite so far on this side as for our 
earlier set of parameters. 



Alternative theory of Grigorov for superconductivity 

Grigorov[Hll [SO] has a completely different type of theory for superconducting 
channels involving pairs of oppositely moving charged density waves along each 
superpolaron string. In his theory he appears to think that the superconduc- 
tivity will occur even for isolated superpolaron strings. If that were the case we 
would have to look for a different explanation for the difference between normal 
and superconducting channels from that used here, and for the difference in 
resistances of the highest-resistance superconducting channels for the two film 
thicknesses. 



Remarks on use of one- dimensional Fermi wave vector 

Our model does not involve a truly one-dimensional system. However, for our 
fitted values of parameters, the transverse wave vector for a given longitudi- 
nal wave number Ul reaches its maximum possible value ir/ax for Ul up to a 
value equal to a large fraction of the calculated 1-D Fermi wave vector for a 
ferromagnetic system of kp,i-D = ttcf = 2ttcb, where cp and Cb are the linear 
concentrations of fermions and bosons. Bearing this in mind, the longitudinal 
Fermi wave vector kpL can be estimated as follows. First we approximate the 
square projection of the Brillouin zone on the rry-plane by a circle of radius k c 
with the same area, i.e. of radius given by 7rfc 2 = (2tt / 1 clt) 2 , or 

k c = 2TT 1/2 /a T . (58) 

Then the maximum value of k m of k at which the transverse wave vector equals 
tt/ot is given by 

k m = (k 2 FL ~pk 2 c ) 1 / 2 . (59) 

If the longitudinal wave vector is kFL, then the volume of the Brillouin 
zone is the same as that for an area irk 2 = 4(7t/q,t) 2 and length 2[(l/4)fcpL + 
(3/4) k m ]. The factor of 3/4 for the weighting of k m occurs because of the 
quadratic dependence of the relevant cross-sectional area on (fc^L — k) up to 
k = k m - Since the volume of the Brillouin zone Vbz — Stt 3 hf — 167r 3 nB, 
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where np and ns are the concentrations of fermions and bosons, we deduce an 
equation for UfLi 

2[(l/4)fc FL + (3/4)(4 L - 4 P 7r/4) 1/2 ](2VaT) 2 = 16^ 3 (c B a T )/4- (60) 

Putting in parameters from our model into these equations, we find that 
kF,i—D — 2Ti(cB<iT)/aT = 2.62 nm _1 and Ufl = 2.95 nm _1 , i.e. the longi- 
tudinal Fermi wave vector is 13% larger than the 1-D value. This is an error 
comparable with or smaller than most of our estimated errors in calculated 
relations between various quantities in our theory. Note that, even if the lon- 
gitudinal and transverse masses were equal, the value of the longitudinal Fermi 
wave vector would be only 37% larger than the 1-D value for the same values of 
other parameters used, and so it is the cut off in transverse wave vectors rather 
than the smaller transverse mass which is the main contributor to the fact that 
the 1-D Fermi wave vector is a fair approximation. 

Effect of replacing the probable tight-binding type of transverse dispersion by one 
with a constant mass 

For a given mass at the bottom of a band, the single-particle bandwidth in the 
x or y direction is larger by a factor of 7r 2 /4 ps 2.5 for a band with a constant 
bare mass than for a tight-binding band. For a pair band made up from single- 
particle states with a tight-binding dispersion, it is less obvious what happens 
to the bandwidth. If we form pair states from two single-particle states with 
equal wave vectors (which may be a fair approximation for weakly bound pairs) , 
but do not consider pair states with transverse wave vector greater than tt/cit 
because of instabilities (cf. [SSI [70]), then the width in the x and y directions 
due to a tight-binding model will only be decreased by a factor (it 2 /8) « 1.23. 
However, whatever happens to the pair bandwidth, if the condensation tem- 
perature is such that ksTs is large compared with the transverse bandwidth, 
we expect the bandwidth is what is most important, whereas for the opposite 
inequality the states near the bottom of the transverse band are what mat- 
ters. For our parameters the transverse bandwidth in the x and y directions 
is p(h 2 /2M L )(it/a T ) 2 = 0.0293 eV, and in the [110] direction twice this. Thus 
we expect that with a tight-binding model the transverse mass at the bottom 
of the band to fit the data would be somewhat smaller that determined by our 
data fitting, but probably not much smaller. 

Value of /ib at T c compared with (iso at T = 0. 

The value of the sound velocity v s at T = in the Bogoliubov theory is 

v s = (mso/Mb) 1 / 2 , (61) 

where hbo is the bosonic chemical potential at T — and Mb is the boson mass 
(see e.g. Ref. [84]). For the sound velocity along the channel direction we use 
the longitudinal boson mass Ml, and we also note that, in the BEC regime[8"2"] 

v s = (0.165 ± 0.085)u F (62) 
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(see section 10), where the Fermi velocity in the longitudinal direction is given 

by 

v FL = hk FL /m FL w 2hk FL /M L . (63) 
We deduce from Eqs. (61) to (63), with Mb = Ml and v F — v F l that 

fi B0 « 0.22(fi 2 k FL /2M L ). (64) 

Using the 1-D value for k F L and putting in values from section 13 for Ml, 
cbcit — 0.259(o,t/p) and a F , we deduce that fiso ~ 0.024 eV. Hence, with our 
value for /is at T c of 0.00127 eV, we find pb/^bo ~ 0.05, with quite a large 
uncertainty because of the uncertainties in Eq. (62). Interpreting the above 
ratio as the fraction of fluctuating condensate at T — T c , we deduce that this 
fraction is about 5%. We do not have a theory for this fraction at present, and 
so cannot say whether its value is what should be expected. 

Coincidence of maximum T c at low currents being only slightly above room tem- 
perature 

Within our theory as presented, the only way to increase the maximum T c would 
be to increase ATi to considerably higher values than the value of 15 K assumed, 
and adjusting other parameters accordingly. This would require ignoring our 
estimated order-of-magnitude value of ATi based on comparison with analysis 
of Zhao for carbon nanotubes assumed to be in the BCS-regime. 

Shortage of experimental input 

Although we have managed to determine twelve parameters in our theory ap- 
proximately, most of the information needed to determine them has come from 
theory rather than from experiment. The only use of experimental data has 
been (i) the use of the histograms in Fig. 1 to determine the probable value of 
the lowest n F for two film thicknesses for which the superconducting T c reaches 
room temperature, (ii) the very approximate estimate of carrier concentration 
in the channels used in section 5, and (iii) the even more approximate estimate 
of the critical current density given in section 10. There is further informa- 
tion from magnetic properties which could be used in principle. However, we 
have had to abandon a theory for the metamagnetic transition given in Ref. 
|52j because recent work [53] on bipolaron dispersion at intermediate coupling in 
one dimension does not support the theory of great enhancements of electron- 
electron attractions at high drift velocities[51 on which the theory was based. 
The large diamagnetism in some samples has been interpreted in terms of in- 
duced currents in superconducting channels forming closed loops [31], but this 
interpretation does not help to determine any of our parameters, except to give 
an alternative estimate of channel diameters to those given in [5] . In Ref. [3T] we 
estimated that diameters for channels which formed closed loops in one sample 
must be greater than 1.52 /im. 
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Some suggested experiments 

The greatest need in this field is for more experimental groups to try to repro- 
duce all experimental results previously reported mainly by one group, at least 
as far as claims of room-temperature superconductivity are concerned, and to 
perform various new experiments. Nine suggestions for experiments were given 
in section 3 of Ref. 52 . Two of these were mainly to test a model for metam- 
agnetism given there which we no longer believe, and another suggestion may 
not be practicable, but most of the others seem to be worthwhile. I mention 
these here more briefly than in Ref. [52 , and also make a couple of new sug- 
gestions. Those previously suggested which still seem to be worthwhile are: 1. 
Use micro-Hall probes near the sample surface to measure the time and spatial 
dependence of magnetic fields near the surface. One property which would be 
determined from such measurements is the direction of the magnetic moment in 
the channels in metamagnetic samples; 2. Perform magnetisation measurements 
for fields parallel to sample surfaces. One might expect a decrease of magneti- 
sation with time in such fields as opposed to an increase seen in fields perpen- 
dicular to surfaces, attributed to gradual alignment of magnetism in channels 
in the direction of the film thickness [19]; 3. Attempt to study the channels by 
optical microscopy at oblique incidence in order to have a component of electric 
field parallel to the channel directions, and with various wavelengths to try to 
find the size of any energy gap in the electron spectrum; 4. Perform electri- 
cal measurements in magnetic fields higher than those required to produce a 
metamagnetic transition to see whether the high conductivity remains in this 
region; 5. Perform magnetic measurements below the glass transition temper- 
ature to see if time dependent effects observed above this temperature [19] in 
magnetic fields no longer occur at low T; 6. Perform electrical measurements 
with films under tension between surfaces to see whether this enables one to 
find conducting channels parallel to film surfaces. 

Another suggestion which would help to test our model, which involves resis- 
tance mainly due to phase slips in the smaller-resistance mediumly conducting 
channels in Fig. 1, would be first to determine histograms of resistance for micro- 
contacts on a suitable film at room temperature, and then to look at the temper- 
ature dependence of resistance at points with resistance near the high-resistance 
end of the range for superconducting channels and at the low-resistance end of 
the range of the mediumly resistant points. If our interpretation, which in- 
volves some resistance due to phase slips, is correct, then the lower-resistance 
mediumly conducting points would become highly conducting at lower temper- 
atures, while points at the higher-resistance end of the low-resistance points 
would change to medium resistances at temperatures somewhere between 295 
K and (295-ftf + AT) , where AT corresponds to the change of transition tem- 
perature due to phase slips in the highest-resistance superconducting channels 
for whatever film thickness is being studied. Quantitative analysis of such ex- 
periments would enable AT and T c to be determined as a function of nr, and 
so would remove our requirement to make guesses at AT. It would also be 
useful to determine the relation between size of channels and their resistance 
(in two-probe measurements) if the size of the channels can be determined e.g. 
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by electron microscopy. Although the experiments suggested in this paragraph 
would be quite time comsuming, we think they would go a long way towards 
establishing superconductivity at room temperature, and whether our model 
for it is realistic. High resolution electron microsopy could possibly detect the 
nm-size subsystems of conducting channels used in our model. 

New theory needed 

Two types of theory which would be useful to develop are: (i) An extension of 
the Ambegoakar-Langer-McCumber-Halperin theory to help determine amounts 
of phase slips for Bose-gas superconductors; (ii) Calculations of transverse ex- 
change energies and any magnetic anisotropy in our model. The first of these 
types of theory is likely to be difficult. To calculate the transverse exchange en- 
ergy using Grigorov's theory [48^ and our parameters should be straightforward, 
but calculations of magnetic anisotropy and sizes of domains do not appear to 
be simple. 

15. Conclusions 

We have used a model of Bose condensation of bosons in arrays of nanofila- 
ments to interpret some results on room-temperature superconductivity in nar- 
row channels through films of oxidised atactic polypropylene. The model makes 
use of a theory of Grigorov and coworkers involving nanofilaments composed 
of strings of an unusual type of polaron produced by interaction with rotat- 
able dipolar groups in an elastomer, and then assuming that the polarons in 
the strings bind to form bipolarons due to mediation of high-frequency bosons. 
These bosons are assumed to be plasmons for most of our calculations, although 
a fair fit to experimental and theoretical constraints may also be obtained assum- 
ing the bosons are high-energy phonons with energy 0.38 eV. There are several 
differences from a model previously used[3], especially (i) a different method of 
estimating the numbers of nanofilaments in the highest-resistance superconduct- 
ing channels, and (ii) a smaller slope of the initial linear term in the dispersion 
of the bosons than previously used. This slope is based on a Bogoliubov form 
of dispersion for bosons at T c with a slope determined by requiring that the 
initial change from the condensation temperature for an ideal Bose gas as the 
boson concentration is increased is about the same as determined by accurate 
published calculations of the effects of interactions on Tb of a Bose gas. There 
are twelve parameters in our model. For one parameter related to changes in 
T c due to phase slips we can only determine a lower limit, and we choose a 
value slightly higher than this limit. The other parameters can then be deter- 
mined approximately by use of twelve relations between parameters of varying 
degrees of estimated accuracy, four of these involving experimental input and 
eight based on theoretical relations between parameters. Some suggestions for 
further experiments are made, and also for theory which it would be worthwhile 
to develop. 
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